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An appropriate extension of the effective potential theory is presented that permits 
the approximate calculation of the dynamical correlation functions for quantum 
systems. These are obtained by evaluating the generating functionals of the Mori 
products of quantities related to the relaxation functions in the (PQSCHA) pure 
self consistent harmonic approximation. 



1 Dynamic Correlations and Spectral Shapes 

The thermodynamicELof quantum systems has been widely studied by the effec- 
e potential theoryEJ: equilibrium properties have been accurately determined 
but the investigation of dynamical quantities is much more difficult . It 
is our purpose to give some insights on how the dynamical problem can be 
approached in terms of the Kubo_.relaxation functions, that naturally appear 
in the framework of Mori theorya. These functions are obtained by suitably 
defined scalar products, the Mori products 

R A , B {t) = (A\B(t)) = / du(A(0)B(t + iu)), (1) 



where the observables A and B are taken such that (A) = (B) = 0, and 
braces denote the thermodynamic average. The Laplace transform of Eq.(Q) 
and in particular, the self-relaxation function So(i) = (Fo\Fo)~ x (Fo(t)\Fo) 
of an hermitian operator Fo(t) can be Laplace transformed and expanded in 
a continued fraction, namely ^j(z) = (z + Sj + i'Ej + i(z))^ 1 , where Sj + i = 
(Fj iFj)^ 1 (Fj + i \Fj + i) and Fj denotes the so-called j-th fluctuating force. The 
quantities (Fj\Fj), can be related to a combination of the first 2(j + 1) mo- 
ments of the time series-expansion of So(f). While (Fj\Fj) with j ^ can be 
expressed in terms of static correlations of time derivatives of Fq , the quantity 
(Fo\Fq) requires the direct evaluation of the Mori product. 
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Indeed, experiments measure the spectral shape, related to S (z) by the 
"detailed balance" principle: 



= (F \F Q ) 



- ljft(5 (z = ia;)). 
1 — e p u 7r 



(2) 



Therefore, <S(w) can be approached from the knowledge of the static quantities 
5j up to a sufficiently large number j = Jtj, Sj(t)uLrul. 



2 Imaginary-time Ordered Products 

Here we shall provide a rigourous derivation of imaginary time ordered prod- 
ucts which allows us to control also the evaluation of real time correlators. 

We start from the generating functional in the hamiltonian path-integral 
form: 



Z [L,J] = j>V[x{u)} J V[p(u)} exp[-±S] 



S = 



dul—ip(u)x(u) + Tt(p(u), x(u)) — KL(u)p(u) — KJ(u)x(u) J . (3) 



According to the effective potential method 0H we use a quadratic trial 
action. The effective Hamiltonian is 7~C elf (ri,£) = w(r],£) + ln(/ _1 sinh/) , 
where / = PHlo/2. By defining the two-component vectors p = (f],Q and 
K(u) = (L(u), J(u)), the approximated generating functional can be written 
as: 



Z [K] 



drj d^ 



2i:h 



exp 



fir* 



fin 



;3h 



du t pK{u) + -j du j dv t K(u)<S>(u-v)K(v) .(4) 



In Eq.(|) we have introduced the 2x2 matrix &ke(u — v) with elements 
$11 (u— v) — m 2 bj 2 <&22{u—v) — m 2 u 2 Kf{u—v) and <J>i 2 (u— v) = — <f> 2 i(w— v) — 
Tf(u — v) , where 



A f {u-v) 
T f (u-v) 



h 



cosh( \ui(u — v)\ — f ) 



2mu> sinh / 

ih sinh( \uj(u-v)\ — f ) 
"Y sinh/ 



sinh / 



(u — v ) — 6(v — u) 



(5) 



While Ay is always well defined, the value of T f for u = v is determined only 
when the limit v — u — > ^ is specified, reflecting the commutation relation of 
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x and p at the same time. Moreover we notice that, A/(0) = Af(f3h) = a and 
that Af(u — v) and Tf(u — v) have a vanishing average in [0, (3h\. 

Defining the two-component vectors z = t (p, x) and y = (yi(u), y 2 ( u ))j the 
following general formula can be derived in the low-coupling approximation! 

(6) 

Where, M = Z 1 [0] is the normalizing factor. The Gaussian average of the 
variables y^ {u v ) is defined by the second moments 

=$«(«-«)• ( 7 ) 

3 Discussion. 

The key result i is represented by the expression (^), with the definition (^). 
Indeed, complicated static Mori products, i.e. moments of any order, can 
be evaluated by this last equation. Static correlations can also be obtained 
performing the appropriate limit u — v — > 0. nn 
For instance, when [u — v) — » + , we indeed recover our previous result aa 
for static averages: 

{ AB)=N /M^je-W.feO, (8) 

J ZTTh 

As far as the Mori product is concerned, the well known series expansion 
(2 7 re)- 1 /2 exp{-x 2 /(2e)} = J2n=o (V«0 [e/2) n {x) appears to be an effi- 
cient tool to approximate the static Mori product of general operators, as well 
as their dynamical correlations, when the scales of the quantum fluctuations 
in the system, ruled by % and the natural length scale a, are small. It has to 
be noticed that the averages of A/(u) and Tf(u) in [0,ph] are vanishing. An 
expansion in terms of these quantities has been recently developed 0. 

At the lowest order (A(j>, x)\B(p, x)) reduces to the "classical like" average 
of the product of the Gaussian spreads of the two operators taken at the same 
order: 

(A(p, x)\B(p, x)) = N(3h (^ ((A( v , 0)) ((B( V , £)» (r >' ^ +o(a, h) . 

J 2irn 

n(9) 

This zeroth order coincides with the assumption proposed in an attempt Q to 
approach also the quantum correlators. 
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In order to do this, let us notice that harmonic oscillators evolve by the 
same law both in classical and in quantum dynamics: the differences between 
quantum and classical statistical evolution are due to the thermal occupation 
numbers, that are static quantities. Let us suppose that the system evolve 
with our effective Hamiltonian as found for the thermodynamic behaviour. The 
Weyl-representationB gives us an unified scheme for describing the dynamical 
variables. The coupling constant g rules the quantum deviations from the 
harmonic behaviour and for vanishing g, Eq.(|^) can be assumed to maintain its 
validity at different real times, provided also that a is small enough; the result 
is exact when g — ► 0. For finite values of g, the validity of this scheme involves 
also the amplitudes of the Gaussian fluctuations ruled by the parameter a. 
Therefore, there is the same behaviour found for approaching static correlators 
with the effective HamiltonianB. 

At this level, the averages in time of tha-quantities ((A)) and ((B(t))), 
evolving with the effective Hamiltonian TL Qit can provide an approximation 
for the time-dependent Mori product (A\B(t)). 

This procedure yields a good approximation for times up to the order of 
7i/3, for which the use of the effective potential makes sense in the calculation 
of the static quantities at lowest order, reproducing for instance a correct 
second moment for the displacement-displacement dynamic correlator with a 
well-behaved classical long time decay. 
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